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TITLE OF THE INVENTION 

PROTECTION RELAY 

CROSS-REFERENCE TO RELATED APPLICATIONS 
This application is based upon and claims the 
5 benefit of priority from the prior Japanese Patent 

Application No. 2001-057948, filed March 2, 2001, the 
entire contents of which are incorporated herein by 
reference . 

BACKGROUND OF THE INVENTION 
©10 1. Field of the Invention 
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IR The present invention relates to a protection 

H 

13 relay in which an influence of distortion component of 

, fault current generated by a fault in a power system is 

suppressed. 
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J? 15 2. Description of the Related Art 
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Generally, a protection relay is used to monitor 
a power system. A main technical subject of such a 
protection relay is to reduce an influence of harmonics 
in fault current and fault voltage generated in a fault 
20 in a system contained in a signal inputted from the 

system. 

Particularly because charging capacity of a system 
is increased in underground power transmission cable, 
phase modifying capacitor and the like, the order of 
25 generated harmonics tends to be lowered. 

Thus, if it is intended to secure a desired 
damping amount according to a conventional method for 



damping the harmonic component with a digital filter, 
filter delay time needs to be prolonged, so that relay 
operation time is delayed. 

For the reason, an approximation method, which is 
not affected by harmonic theoretically even if such 
harmonics is contained has been employed in recent 
years . 

An example of such approximation method will be 
described with reference to FIG. 1. Power transmission 
impedance constant up to a fault point F will be 
considered with reference to FIG. 1. Voltage and 
current of a protection relay installation point A 
are assumed to be v and i when resistance is R and 
inductance is L. If it is assumed that the voltage 
at the faultal point F is zero, differential equation 
of a power transmission line 2 can be expressed in 
an expression (1) . 

V = R- I + L- (di/dt) — (1) 

By calculating a differential item (di/dt) of the 
expression (1) approximately/ detection accuracy can be 
improved even if harmonic is not removed with a filter. 
An example of a concrete method for digital calculation 
actually adapted is shown below. 

V m +V m-1 = R- ( i m +i m-l) +L ' dm-im-l) 

V m-1 +V m-2= R ' <im-l +i m-2> +L ' Um-l-im-^) "'< 2 ) 
When a reactance value X (= coq*L) is calculated 

from the expression (2) , the inductance is expressed in 



expression (3) . 

L m = x m / <°0 = 
( v m + v m-l) ; frm-l + i m-2) " ( v m-l + v m-2) ; (jja + im-l) 
(im ^ im-l) ■ (im-1 + i m-2> " (im-1 ~ i m-2) ' (im + im-l> 

-(3) 

L m /L (true value) is as expressed in the 
expression (4) under conditions of the expressions (5) 
and (6), so that frequency characteristic of X m /X (true 
value) is as indicated with the dotted curve of FIG. 2. 

L m /L(true value) =tan ( co 0 T/2 ) /tan(coT/2) — (4) 

Where, i m =Isin ( co t m ) ^ v m =Vsin ( co t+ 6 ) 

i m -i m -l=2Isin( coT/2) cos ( co t m - coT/2) — (5) 
v m +v m _ 1 =2Vcos ( coT/2) sin( a)t m -coT/2+ 6 ) — (6) 

Therefore, frequency characteristic of X m /X (true 
value) is as indicated with the dotted curve of FIG. 2. 

In FIG. 2, its abscissa axis indicates frequency 
(order) and its ordinate axis indicates a reactance 
measurement value when basic frequency of system 
electrical quantity is 50 Hz. Further, in FIG. 2, its 
dotted curve indicates a case of sampling at 600 Hz and 
its solid line indicates a case of sampling at 4800 Hz. 

As indicated in FIG. 2, the value of L m /L (true 
value) decreases below 1 as the frequency departs from 
its fundamental wave. FIG. 2 indicates that the value 
of (coT/2) only should be suppressed to substantially 1 
(that is, the sampling period should be set small) when 
this value (L m /L) is near twice or three times the 



fundamental wave. 

Frequency characteristic when the sampling 
frequency is actually multiplied eight times is 
indicated by the solid line of FIG . 2. Qauntitavely, 
the relation between an approximate amount (i m - im-l) 
of the differential item and a differentiated amount 
(v m +v m _i) is indicated by the expressions (7) and (8) . 
Therefore, if the sampling frequency is raised (the 
period is decreased) , approximation accuracy of the 
differential item can be raised. 

sin(coT/2) = ooT/2 N cos(coT/2)=l 

i m -i m -. 1 =2I • sin ( coT/2) - cos ( oo t m - co T/2 ) 

=21- co T/2 -cos ( co t m - co T/2) (7) 

v m +v m _ 1 =2V-cos(coT/2) • sin ( co t m - <o T/2 + 0 ) 

=2V'sin(cot m -coT/2+0 ) —(8) 

However, the value of the expression (7) is a very 
small value with respect to an amplitude value I. 
Therefore, a relative value of noise (quantization 
error generated at the time of A/D conversion, white 
noise generated in an analog circuit) contained in 
sampling data (i m , i m -]_) is increased thereby disabling 
practical use of this method. 

For example, when the sampling period T is 
T = 1/4800 sec, coo = 271-50 Hz, the second item 
s/(coOT/2) on the right side of the expression (9) is 
amplified to about 30 times. An s in the expression 
(9) is noise error. 



(i m -i m -l)/(<»0T/2) 

=21- (o)/co 0 ) -cos(cot m -coT/2) + * /(co 0 T/2} --(9) 
BRIEF SUMMARY OF THE INVENTION 

Accordingly, an object of the present invention is 
to provide an improved protection relay* 

A protection relay for determining whether or not 
a faultal point of a power system exists in a predeter- 
mined range, according to a first aspect of the present 
invention is characterized by comprising: filter means 
for inputting sampling data of a voltage and a current 
in the power system to a digital filter having a 
predetermined transfer function and outputting a first 
voltage data and a first current data, and a second 
voltage data and a second current data normal to 
the first voltage data and the first current data, 
respectively; calculation means for calculating 
a predetermined measurement value based on the first 
voltage data, the first current data, the second 
voltage data and the second current data at a first 
time and the first voltage data, the first current 
data, the second voltage data and the second current 
data at a second time different from the first time; 
and operation decision means for performing 
an operation decision based on the predetermined 
measurement value obtained by the calculation means* 

A protection relay for determining whether or 
not a faultal point of a power system exists in 



a predetermined range, according to a second aspect of 
the present invention is characterized by comprising: 
filter means in which sampling data of voltage and 
current in the power system is inputted to a predete- 
rmined transfer function so as to output first voltage 
data and first current data and second voltage data and 
second current data normal to the first voltage data 
and the first current data, respectively; polarized 
voltage value calculation means for inputting the first 
and second voltage data and the first and second 
current data so as to calculate a third voltage normal 
to the first voltage; and operation decision means for 
performing an operation decision based on the third 
voltage . 

A protection relay for determining whether or not 
a faultal point of power system exists in a predeter- 
mined range, according to a third aspect of the present 
invention is characterized by comprising: first filter 
means for inputting sampling data v m and i m of voltage 
v and current i in the power system to a digital 
filter having transfer function f (Z) • (1+k- Z _1 +Z- 2 ) 
(Z indicates a Z conversion operator) so as to output 
voltage data v sm and current data i sm ; second filter 
means in which the sampling data v m , i m are inputted to 
a digital filter having transfer function f(Z)-(l-Z~ 2 ) 
(Z indicates a Z conversion operator) so as to output 
voltage data vj m and current data ij m normal to the 



voltage data v sm and current data i sm ; charging current 
compensation calculation means for calculating quantity 
of electricity defined in i S m" c ' v jm b Y the current 
data i sm , the voltage data Vj m , and a setting value 
C s at time t m ; transmission and reception means 
for transmitting output of the charging current 
compensation calculation means to an opposite terminal 
and when quantity of electricity at the opposite 
terminal is assumed to be B, receiving quantity of 
electricity defined by (i S m~ c " v jm) B at the opposite 
terminal; and operation decision means for performing 
an operation decision based on outputs from the 
charging current compensation calculation means and 
the transmission/reception means according to the 
following expression: 

|| (i S itr Cs - v jm> + (ism- Cs - v jm) B II = 

ka- { II i S m-Cs-Vj m || + || (i sm -Cs • Vj m ) B || } + kb 
where, ||am|| represents a quantity parallel to amplitude 
of AC quantity of electricity "a" at time t m ; 

ka represents a proportion restricting 
coefficient; and 

kb represents minimum sensitivity current. 

Additional objects and advantages of the invention 
will be set forth in the description which follows, and 
in part will be obvious from the description, or may 
be learned by practice of the invention. The objects 
and advantages of the invention may be realized and 



obtained by means of the instrumentalities and 
combinations particularly pointed out hereinafter . 
BRIEF DESCRIPTION OF THE SEVERAL VIEWS OF THE DRAWING 

The accompanying drawings, which are incorporated 
in and constitute a part of the specification, 
illustrate embodiments of the invention, and together 
with the general description given above and the 
detailed description of the embodiments given below, 
serve to explain the principles of the invention. 

FIG. 1 is a power system diagram, which is 
an object of the present invention; 

FIG. 2 is a diagram showing frequency 
characteristic of a reactance measurement value which 
is an object of the present invention; 

FIG. 3 is a diagram showing frequency character- 
istic of a reactance measurement value in a digital 
filter structure, which is an object of the present 
invention; 

FIG. 4 is a block diagram of a protection relay 
according to a first embodiment of the present 
invention; 

FIG. 5 is a block diagram of a protection relay 
according to a first modification of the first 
embodiment; 

FIG. 6 is a diagram showing reactance 
characteristic on an impedance plane; 

FIG. 7 is a block diagram of a protection relay 



according to a second modification of the first 
embodiment; 

FIG. 8 is a block diagram of a protection relay 
according to a third modification of the first 
embodiment; 

FIG. 9 is a diagram showing Ohm characteristic on 
an impedance plane: 

FIG. 10 is a block diagram of a protection relay 
according to a second embodiment of the present 
invention; 

FIG. 11 is a diagram showing offset mho 
characteristic on an impedance plane; 

FIG. 12 is a block diagram of a protection relay 
according to a third embodiment of the present 
invention; 

FIG. 13 is a diagram showing the relation between 
maximum sensitive angle of mho characteristic and 
permanent impedance ; 

FIG. 14 is a diagram showing mho characteristic 
expressed with reference to current; 

FIG. 15 is a block diagram of a protection relay 
according to a first modification of the third 
embodiment; 

FIG. 16 is a block diagram of a protection relay 
according to a second modification of the third 
embodiment; 

FIG. 17 is a block diagram of a protection relay 



according to a fourth embodiment of the present 
invention; and 

FIG. 18 is a diagram for explaining a telegraph 
equation of a power transmission line. 

DETAILED DESCRIPTION OF THE INVENTION 

Prior to description of respective embodiments of 
the present invention, the principle of the present 
invention will be described. If sampling data of 
current i = I -sin (cot) is passed through transfer 
function f(Z) = — +Z~ n ) , current i f sm at 

time t m indicated in the expression (10) is obtained. 
1 1 sm =I • (sin ( co t m ) + sin ( co t m - coT) + ---+sin ( co t m -nco T) ) 

=1- (sin( (n+1) coT/2) /sin(coT/2) ) • sin ( co t m -n co T/2) 

- (10) 

If this current i f sm is passed through transfer 
function (1+k- Z -1 +Z" 2 ) , current i sm at time t m 
indicated in the expression (11) is obtained. 

^-sm"^ ' sm + k* i T sm-l + i 1 sm-2 

=1- (k+2cos(coT) ) • (sin( (n+1) coT/2)/ 

sin (co T/2) ) -sin(cot m - (n+2) coT/2) — (11) 
Likewise, current i = I • sin (ot) is passed through 
the transfer function f (Z) • (1-Z~ 2 ) , current ij m at time 
t m is obtained as indicated in the expression (12) . 
ij m =2I - cos ( co T/2) sin( (n+1) coT/2) - cos ( o>t m - (n+2) coT/2) 

-(12) 

Voltage can be obtained in the same manner as 
current and if it is assumed that the voltages are v sm 



and Vj m , they are in such a relation that they are 
normal to each other. 

Thus, if as the fundamental wave, its value which 
is as near 1 as possible is selected as long as the 
quantity for determining the size of the right side 
of the expression (12) is not | sin ( (n+1) ©T/2) | « 1, 
performance of frequency characteristic can be assured 
while suppressing amplification of noise error. 

Therefore, if ij m/ i S iri' and v sm at time t m and 
ijm-l' ism-l, and v sm-l at time t m _]_ are substituted 
for : 

v sm ' ism-p + ism ' v sm-p 

A m _ ; ; ; •••(13) 

1 jm " 1 sm™p + 1 jm-p * 1 sm 

where m and p represent sampling time series, 

voltage v sm is provided in the form of: 

v sm = V- (k+2cos ( <oT) ) • (sin( (n+1) coT/2) / 

sin (co T/2) •sin(o)t m +0-(n+2) coT/2) •••(14) 

As a result, a reactance value X m is indicated in the 

expression (15) : 

V(k + 2 cos (©T) ) 

Xm = — — — • sin (9) -(15) 

I sin (coT) 

where 9 represents a leading angle of Voltage to 
current . 

Therefore, when n is increased sufficiently in 
the expression (12), an influence of the noise error 
can be reduced. Frequency characteristic of X m 
relative to the fundamental wave can be indicated in 
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a graph shown in FIG. 3. From this graph, it is 
evident that k = 4 provides the most excellent 
frequency characteristic. 

Xm/X (value at fundamental wave) 

(k + 2 cos (coT) ) / (k + 2 cos (co 0 T) ) ^ 

•••(16) 

sin (coT) / sin (co 0 T) 

A protection relay according to a first embodiment 
of the present invention will be described based on 
the above-described principle with reference to FIG. 4. 
FIG. 4 is a block diagram showing schematically the 
structure of the protection relay according to the 
first embodiment. The protection relay of the first 
embodiment suppresses amplification of the noise error 
contained in minute quantity of data even when the 
sampling period is shortened, in order to secure 
performance of frequency characteristic indicated 
in the expression (4) . 

Referring to FIG. 4, a first (digital) filter 1 
comprises a digital filter (Z indicates a Z conversion 
operator) which outputs sampling data v m and i m about 
voltage v and current i of a power system as voltage 
data v sm and current data i sm through transfer function 
f (Z) • ( 1 + k- Z"' 1 +Z"" 2 ) and extracts predetermined frequency 
components of voltage and current of a power system 
(not shown), which is a protecting object. 

A second (digital) filter 2 comprises a digital 
filter (Z indicates a Z conversion operator) which 



outputs sampling data v m and i m as voltage data Vj m and 
ij m through transfer function f(Z)-(l-Z~ 2 ) and extracts 
voltage and current normal to the first filter 1 in all 
frequency components. 

A reactance value calculator 3 calculates 
a reactance value X m from voltage data v sm , current 
data i sm/ voltage data Vj m and voltage data ij m at time 
t m and voltage data i S m-p' current data i sm -p/ voltage 
data vj m _p and current data ij m -p at time t m _p. 

An operation decision section 4 compares the 
reactance value X m obtained by the reactance calculator 
3 with a preliminarily set integer (whether or not 
X m < X s is determined with X s ) . The preliminarily 
set integer value is called setting value in this 
specification. 

With such a configuration, the first filter having 
the transfer function f (Z) • (1+k* Z~ 1 +Z"" 2 ) receives input 
of sampling data v m and i m about voltage v and current 
i of the power system and outputs voltage data v sm and 
current data i sm . 

Further, the second filter having the transfer 
function f(Z)-(l-Z~ 2 ) receives input of sampling data 
v m and i m about the voltage v and current i and outputs 
voltage data Vj m and current data ij m . 

The reactance value calculator 3 calculates 
a reactance value X m from the voltage data v sm and 
current data i sm obtained by the first filter 1 and 



voltage data Vj m and current data ij m obtained by the 
second filter 2 at time t m , and voltage data v SItl _p and 
current data i S m-p obtained by the first filter 1 and 
voltage data Vj m _p and current data ij m -p obtained by 
the second filter 2 at time t m _p according to the 
expression (17) : 

_ - v sm ' ^srn-p + ism ' v sm-p 
X m „ ■ ■ ■ ■ •••(!/) 

~ x jm * 1 sm-p + 1 jm~p * 1 sm 
where, m and p represent sampling time 

series . 

Further, the operation decision section 4 decides 
whether or not X m < X s is established from the 
reactance value X m and setting value X s obtained by the 
reactance value calculator 3 and when that relation is 
established, it decides it is active and otherwise, 
it decides it is inactive. 

The transfer functions of the first and second 
filters 1 and 2 are expressed in f (Z) • ( 1 + k • Z"" 1 + Z" 2 ) 
and f (Z) • (1-Z~ 2 ) when the Z conversion operator is 
employed. As already described before, the outputs of 
the first and second filters 1 and 2 are normal to each 
other. There is a relation that the output of the 
first filter 1 is delayed at 90° with respect to the 
output of the second filter 2. 

It is permissible to define these transfer 
functions with transfer function 1: f(Z), 
transfer function 2: ( 1 + k • Z _1 +Z~ 2 ) and transfer 



function 3: 1-Z~ Z and form the first filter 1 and 
the second filter 2 by combining those. 

That is, the same output can be obtained by 
passing the input voltage and current through a digital 
filter of the transfer function 1 and then passing its 
output through digital filters of the transfer 
functions 2 and 3. 

The reactance value calculator 3 of FIG. 4 
calculates a reactance value from a protection relay 
installation point in a power transmission line of 
FIG. 1 up to a faultal point according to the 
expression (17) . Assuming that the input voltage 
and current are i = I -sin (cot), and v = V- sin (cDt+0) , 
the expression (17) can be expressed in the 
expression (15) . 

The operation decision section 4 corrects the 
decision expression using the reactance value X m , 
setting value X s and preliminarily set fundamental wave 
constant sin (©o T ) / ( k+2 • cos (o>o T ) ) with respect to the 
reactance value X m calculated by the reactance value 
calculator 3 as indicated in the expression (18) so as 
to decide whether or not it is active. Such operation 
decision is carried out plural times ordinarily. 

X m <X s / (sin (o>o'T) / (k+2 - cos (co 0 -T) ) —(18) 

Consequently, amplification of error due to 
approximation of the differential item is suppressed, 
so that L m /L (true value) becomes a characteristic 



which is as near 1 as possible in a wide frequency 
band. Therefore, it is possible to block an influence 
even when a harmonic component occurs in fault voltage 
and fault current of the power system thereby improving 
reliability. 

According to the above description, the reactance 
value calculator 3 calculates the reactance value using 
outputs of the first filter 1 and the second filter 2 
and the operation decision section 4 decides whether or 
not it is active according to the condition expression 
(18) based on that reactance value. 

The present invention is not restricted to this 
example, but it is permissible for an operating 
restarint quantity calculator 5 to calculate a m and b m 
according to the expression (19) based on the outputs 
of the first filter 1 and the second filter 2 and then 
for an operation decision section 6 to decide the 
operation according to the expression (20) to be 
described later using that result. 

a m = ~ v sm* 1 sm-p +v sm-p * 1 sm 

km^-ijirT i-sm-p+i- jm-p * ism ,# " (19) 

That is, when current i = I-sin(<ot) and voltage 
v = V- sin (GDt+0) are inputted to the first filter 1 of 
the transfer function f (Z) • (1-k- Z -1 +Z~ 2 ) and the second 
filter 2 of the transfer function f ( Z ) • ( 1-Z~ 2 ) , the 
expression (19) is converted to the expression (20), 
providing that f(Z) = (1 + Z" 1 + Z'" 2 + + Z~ n ) . 



a m =IV{ (k+2cos ( o)T) ) 2 - 

(sin{ (n+1) coT/2) /sin ( coT/2) ) 2 } • sin ( 0 ) - sin(pcoT) 
b m =2I 2 {cos ( o>T/2) • sin ( (n+1) co T/2) • 

(k+2cos(coT) ) • (sin( (n+1) <o T/2 ) /sin ( co T/2 ) ) } 
•sin(pcoT) •••(20) 
The operation decision section 6 corrects 
sensitivity constant Kq and reactance setting value 
X s to X s <r- X s / (sin (oqT) / (k+2 • cos (0q t ) ) using a m and 
b m outputted in this way and decides the operation 
according to a decision expression of the 
expression (21) . 

Although this is a different method from 
the above-described method, it has the reactance 
characteristic as shown in FIG . 6 like that example. 

Consequently, amplification of error due to 
approximation of the differential item is suppressed, 
so that L m /L (true value) becomes a characteristic 
which is as near 1 as possible in a wide frequency 
band. Therefore, it is possible to block an influence 
even when a harmonic component occurs in fault voltage 
and fault current of the power system thereby improving 
reliability. 

Further, it is permissible to use an Ohm value 
calculator 7 defined by expression (22) as shown in 
FIG. 7 as such calculation means or an operation 
decision section 8 defined by the expression (25) , 



which will be described later. 

R _ ~ jjm ; v sm-p + v sm ; ijm-p --(22) 
~ i jm * ism-p + ijm-p * ^sm 

With this configuration, when current 

i = I -sin (cot) and voltage v = V-sin(ot+9) are inputted 

into the first filter 1 of the transfer function 

f (Z) • (l + k-Z-i + Z" 2 ) and the second filter 2 of the 

transfer function f (Z) . ( 1-Z~ 2 ) , the expression (23) 

is established, providing that f(Z) = (1+Z~ 1 -HZ~ 2 + ... 

+Z~ n ) . Further, the Ohm value R m is converted to the 

expression (24) because of the relation of the 

expression (20) . 

"Mm* v sm-p +:L jm-p " v sm 

= 2IVcos ( coT/2) • sin ( (n+1) coT/2) • (k+2cos ( ooT) ) • 
(sin( (n+1) coT/2) /sin ( coT/2) ) -cos( 6 ) -sin(pcoT) 

-(23) 

R m =(V/I) -cos ( 6 ) — (24) 

The operation decision section 8 decides whether 
or not a decision expression for the expression (25) is 
established based on the Ohm value R m and the setting 
value R s calculated from the expression (24) and when 
that expression is established, it decides it is 
active . 

R m <R s --(25) 
Consequently, amplification of error due to 
approximation of the differential item is suppressed, 
so that L m /L (true value) becomes a characteristic 



which is as near 1 as possible in a wide frequency 
band. Therefore, it is possible to block an influence 
even when a harmonic component occurs in fault voltage 
and fault current of the power system thereby improving 
reliability. 

Further, it is permissible to use an operation 
restricting quantity calculator 9 defined by the 
expression (26) shown in FIG. 8 or an operation 
decision section 10 defined by the expression (28) 
which will be described later, as such calculation 
means . 

c m =_1 jm' v sm-p +v sm # 1 jm-p 

b m =-ij m - i sm _p+i j m -p* i SItl ••• (2 6) 

That is, when current i = I* sin (cot) and voltage 
v = V*sin(cot+0) are inputted to the first filter 1 of 
the transfer function f (Z) * (1-k- Z~!+Z~ 2 ) and the second 
filter 2 of the transfer function f ( Z ) • ( 1-Z~ 2 ) , the 
expression (26) is converted to the expression (27), 
providing that f(Z) = (1+Z~ 1 + Z _2 + + Z~ n ) . 

m = _i jm' v sm-p +i jm-p" v sm 

= 2IVcos ( coT/2) -sin( (n+1) coT/2) • (k+2cos(a)T) ) • 

(sin ( (n+1) coT/2) /sin ( coT/2) ) • cos ( 8 ) • sin (pco T) 

-•(27) 

Then, the operation decision section 10 decides 
whether or not it is active based on the decision 
expression (28) from c m , b m , Ohm setting value R s 
and sensitivity constant - 



Rg-Cm^k! — (28) 

Such protection relay has Ohm characteristic 
shown in FIG. 9 like the configuration described with 
reference to FIG. 7 before and is different from it 
only in its realization method. 

Consequently, amplification of error due to 
approximation of the differential item is suppressed, 
so that L m /L (true value) becomes a characteristic 
which is as near 1 as possible in a wide frequency 
band. Therefore, it is possible to block an influence 
even when a harmonic component occurs in fault voltage 
and fault current of the power system thereby improving 
reliability. 

Next, a second embodiment of the present invention 
will be described with reference to the accompanying 
drawings. A description of the same configuration as 
the above-described embodiment is omitted appropriately 
while like reference numerals are attached to the same 
components in the Figures. 

According to the first embodiment, for example, 
FIG. 4 indicates calculation of reactance and FIGS. 5, 
8 indicate calculation of operating restarint quantity, 
and FIG. 7 indicates calculation of Ohm value and 
decision of operation based thereon. 

This embodiment contains the Ohm value calculator 
7 and reactance value calculator 3 as shown in FIG. 10 
and then, the operation decision section 11 decides 



the operation using these outputs. 

That is, the voltage v sm , current i sm and voltage 
Vj m , current ij m outputted from the first filter 1 and 
the second filter 2 are inputted to the Ohm value 
calculator 7 for calculating Ohm value and the 
reactance value calculator 3 for calculating reactance. 
The Ohm value calculator 7 calculates an Ohm value R m 
according to the expression: 

_ x jm ' v sm-p + v sm ' ^jm-p 
m ~ ■ 7™ : 

^-jm * 1 sm-p + ^-jm-p ' 1 sm 

The reactance value calculator 3 calculates 

reactance X m according to the expression (13) . 

The operation decision section 11 decides whether 

or not it is active according to the expression (29) 

using the Ohm value R m and the reactance value X m . 

This protection relay has offset mho characteristic 

indicated in FIG. 11: 

(R m -R 0 ) ' (RirrV + (X m -X 0 ) • (X m -X F ) < 0 ••• (29) 

where, Rq (Ohm component) represents an offset mho near 

side setting value; 

Xq (reactance component) represents an offset 

mho near side setting value; 

R|r (Ohm component) represents an offset mho far 

side setting value; and 

Xp (reactance component) represents an offset 

mho far side setting value. 

A third embodiment of the present invention will 
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be described with reference to the accompanying 
drawings . Description of the same configuration as 
the above-described embodiments is omitted while like 
reference numerals are attached thereto* 
5 A protection relay of this embodiment comprises a 

polarized voltage value calculator 12 in which voltage 
v sm , current i sm and voltage Vj m , current ij m are 
inputted in order to extract voltage v p i m normal to 
1^ voltage v sm and an operation decision section 13 in 

|% 10 which the aforementioned voltages and currents are 



18 



inputted while voltage Vpi m is inputted from the 



J£ polarized voltage value calculator 12 so as to 

M 

decide whether or not it is active according to 
the expression (30) . 

3 tt 

O 15 v pjm _ p - ( (Rs-i sm +Xs-ij m )-v sm )- 

O v pjnT ( (Rs * i SIU _p+Xs • i j m _p) ""V sm _p) > K2 -(30) 

ru 

In the meantime, R s and X s are setting value 
s of an Ohm component and a reactance component 
respectively, while X s is corrected in the form of 
20 X s <- X s /(sin(©OT)/(k+2-cos(<oOT) ) . 

In the expression (30) , the quantities of 
electricity in the item of (Rs'^sm + x s -ij m ) and item 
of (Rs'^-sm-p + x s*^jm-p) are ^ n such a relation that 
they are advanced only in size of (R s 2 + X s 2 ) 1//2 and 
25 phase of <f>= tan""! (X s /R s ) with respect to current i sm . 

Its example will be shown in FIG, 13. 

If current i = I • sin (cot) and voltage 
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v = V- sin (oot+0) are inputted to the first filter 1 of 
the transfer function f (Z) • (1-k- Z"" 1 + Z"' 2 ) and the 
second filter 2 of the transfer function f (Z) • (1-Z~ 2 ) , 
the expression (31) is established, providing that 
5 f(Z) = (l+Z- 1 +Z- 2 + — +Z~ n ) : 

V p j * {Rs- I s -cos ( Q )+Xs/ (sin( co 0 T) / 

(k+2cos (o) 0 T) ) ) *Ij -sin( 0 ) ) -Vs} -sin(pcoT) >K2 
where: I S =I (k+2cos (coT) ) (sin( (n+1) ooT/2) /sin (co T/2) ) 
^ V s =V(k+2cos(coT) ) (sin( (n+1) coT/2) /sin ( coT/2) ) 

C3 

S 10 I -a =2 1 cos ( co T/2 ) sin( (n+1) coT/2) (31) 

Iff Then, if the expression (31) is expressed with 

S| 

Q Zs (size: (R s 2 + Xs 2 ) 1 / 2 and phase * = tan" 1 (X s /R s ) , 

Ci 

a expression (32) is established. Therefore, the 

C3 

Q; expression (32) indicates operation principle 

Q 

i'7 15 expression of mho characteristic shown in FIG. 14 if 



ru 



its frequency is basic frequency. 

Vpj • { (Zs • Is - cos ( 6 ) cos ( <£ ) + Zs • Ij • sin ( 0 ) sin ( 4> ) / 

(sin ( co qT) / (k+2cos ( a) qT) ) -Vs) • sin (pcoT) 
=Vpj ■ { Zs • I • (cos ( 6 ) cos ( ) (k+2cos ( coT) ) / 
20 sin(coT)+sin( 0 ) sin( 4> ) 

(k+2cos ( co qT ) ) / sin ( ( co qT) ) -Vs} • sin (pcoT) 

V CO — CO Q 

=Vpj • {Zs-I-cos ( 0 - <f> ) -Vs} -sin(pcoT) >K2 — (32) 

The above-described structure indicates a case 
25 where voltage Vpj m normal to voltage v sm in the 

fundamental wave is extracted by the polarized voltage 
value calculator 12. The present invention is not 
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restricted to this example. For example, it is 
permissible to provide with a polarized voltage 
value calculator 14 which regards a voltage before 
predetermined cycles (data before N samplings) of the 
5 voltage vj m normal to voltage V sm as polarized voltage 

as shown in FIG. 15. In this case, the voltage v p j m 
and voltage v sm are indicated in the expression (33) . 
v p j m =2V'sin( (n+1) coT/2) - cos (coT/2) 
t* -cos (cotm+ 0 =NcoT- (n+2) a>T/2) 

C3 

O 10 v sm =V- (k+2cos (a>T) ) -sin( (n+1) coT/2)/ 

(§ 

m sin ( co T/2) -sin( cotm+ 6 - (n+2) coT/2) (33) 

si 

0 (Nco qT=2 % -M, M is an integer) 

ii 

i " It is permissible to provide with a polarized 

jjj voltage value calculator 15 for extracting voltage Vpj_ m 

15 normal to voltage v sm as shown in FIG. 16 and if this 



o 
ru: 



is for detecting a short-circuit, in case of, for 
example, AB phase, extract positive sequence voltage 
relative to the AB phase. A, B, and C indicate each 
phase of three-phase AC quantity of electricity. 

20 For example, positive sequence voltage relative 

to the AB phase can be extracted according to 
Vpsm^) = -3 1 / 2 (vj m (C) - Vj m (0)) + v sm (AB). Here, 
v sm (0) indicates zero phase quantity of electricity. 

Further, if this is for ground fault, the positive 

25 sequence voltage relative to A phase can be calculated 

according to v psm (A) = 3 1 / 2 - (v sm (A) - v sm (0) + Vj m (BC) . 
In addition to the method for extracting voltage 
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vector from two electricity quantities normal to each 
other, it is permissible to apply a method in which 
sampling time series is shifted by 90° or equivalent. 

Consequently, amplification of error due to 
approximation of the differential item is suppressed, 
so that L m /L (true value) becomes a characteristic 
which is as near 1 as possible in a wide frequency 
band. Therefore, it is possible to block an influence 
even when a harmonic component occurs in fault voltage 
and fault current of power system thereby improving 
reliability. 

A fourth embodiment of the present invention will 
be described with reference to FIG. 17. Description 
of the same configuration as the above-described 
embodiments is omitted while like reference numerals 
are attached thereto. 

FIG. 17 is a block diagram of a protection relay 
applied to description of this embodiment. According 
to the fourth embodiment, outputs from the first filter 
1 and the second filter 2 are inputted to a charging 
current compensation calculator 16 and its output is 
inputted to an operation decision section 18 which 
follows operation decision principle of ratio 
differential relay type so as to decide whether or not 
it is active. 

The charging current compensation calculator 
16 corrects a setting value C s in the form of 



C s <r- C s - (k+2cos (©oT) /sin(o>o T > and calculates 
i sm -C s * Vj m using output current i sm from the first 
filter 1 and output voltage Vj m from the second 
filter 2. In the meantime, C s -Vj m indicates current 
compensation generated by charge capacity C s . 

Then, opposing quantity of electricity 
(i sm -C s - Vj m ) B of an opposite terminal is received by 
transmission and reception section 17 and quantity of 
electricity of itself terminal is transmitted to the 
electric power station. Here, B indicates quantity of 
electricity at an opposite terminal. 

After that, the operation decision section 18 
decides the operation based on vectorial sum of current 
which compensates for charge current at itself terminal 
obtained by the charging current compensation 
calculator 16 and current which compensates for charge 
current at a terminal of the opposite terminal, namely, 
scalar sum of amplitude value of differential current 
and currents which compensate for charge current 
at each terminal, according to expression (34) : 
II ( ism-Cs •v jm )+(i sm -Cs-v jm )B || > 

lea- { || ( ism-Cs -v jm ) II + II (i sm -Cs • v jm ) B || }+kb 

-•(34) 

where, ||am|| represents quantity parallel to amplitude of 
AC quantity of electricity at time t m ; 

ka represents a proportion restricting 
coefficient; and 



kb represents minimum sensitivity current. 

Physical meaning of the charging current 
compensation indicated in the expression (34) will be 
described with reference to a power transmission line 
in FIG. 18. A well known telegraph equation is 
expressed at transmission/reception terminal as shown 
in expression (35) . Then, if Taylor expansion 
approximation is executed on differential current i^n 
(t) with the condition of (x ^ 0), the expression (36) 
is obtained: 

Waveform of forward movement: 

i DF (t)=i s (t- t )+e s (t-T )/z+i R (t+T )-e R (t+x )/z 
Waveform of backward movement: 

i DB (t)=i s (t- r )-e s (t-r )/z+i R (t+r )+e R (t+r )/z 
Differential current: i DD (t) = (i DF (t) + i D B (t) )/2 — (35) 
where, suffix S represents a transmitting terminal; 
R represents a receiving terminal; 
Z represents a surge impedance =; and 
x represents propagation time = 1/ v (LC) . 

i DD( t) = i s( t) - (, / z) ^ + (,2 / z) - (,3 / 31z) ■ ■ 

+ iR(t) _ (x/Z )^ + (x2/ Z )^-(,3/31z)^> 

- (36) 

V t /Z=Ec/2(l/2 charge capacity in entire SR 
interval) 

Even if differential current is extracted with 
only current vectorial sum current (is(t) + iR(t)) at 
the transmission/reception terminal, the aforementioned 



charge current component becomes error current, thereby 
inducing a drop of sensitivity of the differential 
relay. Thus, only a fault current component can be 
extracted by compensating for this, 

i DD (t) = i s (t) - (x / z) + i R (t) - <t / z) - (37) 

dt dt 

In the above description, the transfer 
function f(Z) in the first filter 1 and the 
second filter 2 has been described in the form of 
f(Z) = (l+Z-!+Z^ 2 + — +Z" n ) . It is evident that the 
same effect can be obtained even if the transfer 
function is f(Z) = (1+Z~ 2 +Z~ 4 + +Z~ 2n ) . 

In this case, in the expression (10), coT/2 only 
should be replaced for ©T and in that case, the 
expression (38) is established. Likewise, the 
expressions (8) and (12) only should be replaced. 
For convenience of following description, it is 
described that fl = f(Z) = (1 + Z"' 1 + Z" 2 + ••• + Z~ n ) , and 
f2 = f(Z) = (l+Z- 2 +Z- 4 + — +Z~ 2n ) . 

1 ' sm =I ' (sin(a)t m )+sin(cot m -2coT) + - • + sin ( to t m -2nco T) ) 

= I* (sin( (n+1) coT)/sin(o)T) ) • sin ( co t m -n co T) — (38) 
If fl is replaced to f2, the quantity of data in 
the latter can be reduced assuming that the window 
length of data for use is the same and approximation 
error of integration is of the same value. For 
example, if n = 2 is adapted in case of fl, n = 1 may 
be adapted in case of f 2 . 



Because gain in the first filter 1 and the second 
filter 2 using fl is larger than gain in the first 
filter 1 and the second filter 2 using f2, noise error 
can be compressed by that amount. 

For example, when fl is used, gain Gl of n = 2 is 
expressed in the expression (39) and if f2 is used, 
gain G2 of n = 2 is expressed in the expression (40) . 
Gl= | sin(3coT/2)/sin(coT/2) | = | 2cos(<»T)+l | •••(39) 
G2= | sin(2coT)/sin(o)T) | = |2cos(coT) | --(40) 

Such transfer function f (Z) is constructed with 
FIR fileter. The present invention is not restricted 
to this example, but the same integration error 
characteristic can be achieved even when it is 
constructed with recursive digital filter. 

Consequently, amplification of error due to 
approximation of the differential item is suppressed, 
so that L m /L (true value) becomes a characteristic 
which is as near 1 as possible in a wide frequency 
band. Therefore, it is possible to block an influence 
even when a harmonic component occurs in fault voltage 
and fault current of the power system thereby improving 
reliability. 

According to the present invention, as described 
above, even when a harmonic component is overlaid on 
a fault current and a fault voltage generated upon 
a fault of a power system, a predetermined time 
differential equation can be solved approximately by 



passing through predetermined digital filters normal to 
each other in a wide frequency band, thereby achieving 
a high accuracy protection relay. 

Additional advantages and modifications will 
readily occur to those skilled in the art. Therefore, 
the invention in its broader aspects is not limited to 
the specific details and representative embodiments 
shown and described herein. Accordingly, various 
modifications may be made without departing from the 
spirit or scope of the general inventive concept as 
defined by the appended claims and their equivalents. 



